Flux-flow in d-wacve superconductors: Low temperature universality and scaling 
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We demonstrate that superclean ei-wave superconductors with A^t/Ef ^> 1 display a novel type 
of vortex dynamics: At low temperatures, both dissipative and transverse components of the flux- 
flow conductivity are found to approach universal values even in the limit of infinite relaxation 
time. A finite dissipation in the superclean limit is explained in terms of the Landau damping on 
zero- frequency vortex modes which appear due to minigap nodes in the bound-state spectrum in 
the vortex core. In the moderately clean regime the scaling law at low T and low field is obtained. 
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The Landau description of low-temperature thermody- 
namics and low-frequency dynamics of Fermi and Bosc 
liquids is based on the concept of predominating low- 
energy quasiparticles. In the same way, the low- frequency 
dynamics of vortices in Fermi superfluids and supercon- 
ductors appears to be also determined by low-energy ex- 
citations. For Fermi superfluids and superconductors 
with a gap in the spectrum, the relevant low-energy 
quasiparticles are concentrated in vortex cores. The 
anomalous spectral branch of quasiparticles localized in 
the core, which crosses zero as a function of the orbital 
quantum number, was first found by Caroli, de Genncs 
and Matricon jjj for the simplest case of axisymmet- 
ric vortices in s-wave superconductors. Effects of these 
quasiparticles on the vortex dynamics were considered in 
a number of papers (see, for example, (|||). 

The d-wave superconductivity adds new aspects to this 
problem: due to the gap nodes, gapless quasiparticles ex- 
ist also outside vortex cores. The low-energy fermions 
which are far from vortex cores but still are under the in- 
fluence of the vortex superflow velocity (i) introduce non- 
analytical terms in the free energy of the mixed state of 
a d-wave superconductor Q ; (ii) result in a singularity in 
the vortex density of states || and (iii) lead to a peculiar 
scaling law for thermodynamic and kinetic properties of 
d-wave superconductors at low r«T c in small magnetic 
fields B <C H C 2 J|-01. Here we show that the vortex dy- 
namics is also influenced by these far-distant excitations 
bound to vortices. This results in a new scaling law for 
the Ohmic and Hall conductivities in the mixed state of 
a d-wave superconductor. Moreover, a novel effect ap- 
pears: the vortex motion is accompanied by the Landau 
damping on vortex modes giving rise to a finite dissipa- 
tion even in the superclean limit. 

It is known |^|| that, in a d-wave superconductor, 
impurity scattering broadens the gap nodes to the im- 
purity band whose width is 5 ~ T c exp(— ttAqt) in the 
Born approximation. For temperatures below the im- 
purity bandwidth S, the current response function in 



the Meissner state, j = — Q(oj)A(u>) displays unusual 
behavior: The "conductivity" <r(0) defined as Q(lu —> 
0) = n s e 2 /mc — i(uj/c)o~(0), saturates to a scattering- 
independent universal value |],[^] corresponding to the 
Drude conductivity with the effective relaxation time 
1/t ~ T c . In the present Letter, we demonstrate that 
a similar effect exists in the mixed state of d-wave super- 
conductors: Both longitudinal (Ohmic) and transverse 
(Hall) components of the conductivity tensor approach 
finite universal values even in the limit of infinite relax- 
ation time, it temperatures and magnetic fields are low 
enough. This universal behavior of the flux-flow conduc- 
tivity is realized in superclean systems with the mean 
free path I {Ep/T c )^. We trace the finite dissipation 
in the superclean limit back to the resonant absorption 
(Landau damping) on zero- frequency vortex modes which 
appear due to the nodes in the minigap c^o, i-e., in the dis- 
tance between bound-state levels in the vortex core. Note 
that the width of the abovementioned impurity band S is 
vanishingly small for superconductors with an impurity 
scattering rate satisfying the superclean condition. 

We start with the spectrum of relevant low-energy 
fermionic quasiparticles assuming that T <C T c . Since 
our results are not sensitive to the details of the quasi- 
particle spectrum, we take the simplest spectrum, which 
has the tetragonal symmetry and four nodes. In a bulk d- 

wave superconductor, E(p) = \J £,p + A 2 ^ sin 2 (2a). Here 
£ p = e(p) — Ep. We consider an isotropic Fermi sur- 
face around the c axis and assume that the vortex is 
axisymmetric. The magnetic field is along the symme- 
try axis which is the z axis of the coordinate frame 
with the positive direction along the vortex circulation 
z = sgn(e)B/_B; a is the asimuthal angle of p counted 
from one of the 4 nodes in the gap (a = 0) . 

Quasiparticles bound to the vortex. The low-energy 
spectrum of Caroli-de-Gennes-Matricon quasiparticles 
around a vortex contains an anomalous branch. Exci- 
tations on this branch are characterized by two canoni- 
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cally conjugated variables, the angle a and the angular 
momentum Q. In the quasiclassical approximation the 
variables a and Q are commuting, and the spectrum is 
E(Q,a) — —LUo(a)Q. For a <i-wave superconductor close 
to the gap nodes | sin(2a)| <C 1, the distance between the 
energy levels, i.e., the minigap is ffl5| 



w (a) 
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Here v± are the projections of the Fermi momentum 
and velocity on the x — y plane. According to Eq. (Q), 
the minigap u)q as a function of a has nodes of a higher 
order than the gap itself. We shall show that the strong 
singularity in loq (a) results in a nonvanishing dissipation 
even in the superclean limit A^r/Ep 3> 1 when the 
quasiclassical minigap u>o(a) is well resolved. 

Kinetic equation for the distribution function f(t, a, Q) 
for a system of fcrmions characterized by canonically con- 
jugated variables has a conventional form pjj: 



dt 



d£dE 
da dQ 



dEdf_ 
da dQ 



f-fo 
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Here dE/dQ = — u>o(a); fo is the equilibrium Fermi 
distribution; the collision integral is written in a r- 
approximation. A rigourous consideration shows that the 
collision integral for localized excitations does not have 
singularities as a function of a, thus it can be considered 
as a constant on order of the relaxation time r in the 
normal state. 

If the vortex moves with a constant velocity u = 
u(xcos0 + ysin0) with respect to the lattice, the 
Doppler shift of the energy E — > E(Q, a) — pu pro- 
duces the "driving force" (dE/da) — up± s'm(a — </>) 
acting on quasiparticles. The third term in the lhs of 
kinetic equation (||) contains this force multiplied by 
dfo/dQ — —uio(a)(dfo/dE). Introducing the longitudi- 
nal and transverse response to this perturbation and tak- 
ing into account that, for the tetragonal symmetry, the 
response does not depend on the direction of the vortex 
motion with respect to the crystal lattice, one obtains 



c r $0 r / 

bf = -up± — [^f H cos{a 



+ 70 sin(a - </>)] (3) 



where 70(a) and 7jj(a) satisfy two coupled first-order 
differential equations 
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da 



jh - U(a)-f + 1 = 0, 



da 



+ 7o -U{a)j H = 



(4) 



where U(a) — [wo(a)r] . These two real functions 
are combined into one complex function 70 = ImW^(a), 
7# = ReW(a), which obeys the equation 



(dW/da) - (U + i)W + i = 0. 



(5) 



The force acting on the vortex from the environment 
where it moves is the momentum transferred from the 
excitations: 



dp z dQ da dp 
~2tT 2vr ~dt 



(6) 



Using dp/dt = —S7E — [z x p](dE/dQ) we get 



dp z dQ da . „ . dE 



Z X u 



(7o) F u) . 
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Here n is the electron density; {■ ■ -) F is the average over 
the whole Fermi surface. The force F onv should be bal- 
anced by the Lorentz force = irn[v s x z]. The so 
called mutual friction parameters (70) p and ('Jh)p are 
thus coupled to the flux-flow longitudinal and Hall con- 
ductivities. At low temperatures, 

(lo)p = (Ba /n\e\c), (jh) p = {Bcr H /nee). (8) 

If the Fermi surface has electron-like and hole-like pock- 
ets, the conductivities at low temperatures are p| 



CO 



\e\c 
B 



0-H = - 



n e (70) F:G + n h (70} Fjh J , (9) 
f ne{lH) F . e -n h ( lH ) Fh . (10) 

where {■ ■ ■) F e and (...) Fh are the averages over the 
electron-like and hole-like parts of the Fermi surface, re- 
spectively; n e and rih are the numbers of electrons and 
holes. 

Solution of the kinetic equation (|^) is 



\V = <■'" n ' !i (C-ij e -t <Q ' +F (°')l da' 



(11) 



where F(a) — U{a') da'. In the moderately clean 
limit, A^t/Ep <C 1, the potential U(a) is always large, 
and we obtain the local solution as in an s-wave super- 
conductor pi 



70(a) =u! (a)r, 7if(a) = [w (a)r] z 



The conductivities are 

nec A^r 
B E F 



a~0 



In- 



T 



(12) 



(13) 



and o-h/o~o ~ {^to T /BF ) ln(T c /T) which are similar to 
the results for an s-wave superconductor P,p|,^2]| . 

In the superclean limit, Aqt/E f 1, the situation 
is more interesting. The potential U(a) is small al- 
most everywhere except for vicinities of the gap nodes 
where it diverges: U(a) ~ l/[r(<5a) 2 ] at \Sa\ <C l.Here 
r ~ Aqt/E f . However there is a physical infrared cutoff 
drain for the angle below which the divergent potential 
U(a) saturates. For example, if the vortices are separated 
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by a finite distance R ~ ^ H c2 / B with B <C H C 2, i.e., 
R^> the Doppler shift of the energy pfV s ~ T c (£/R) 
(where v s is the vortex-induced superflow velocity) af- 
fects the spectrum Eq. (|l|) when the excitation energy 
becomes comparable with the gap Ao| sin(2a)|. This de- 
termines the cutoff a rn i n = ^JB/H c2 . In addition, a 
finite temperature introduces a cutoff a m i n = T/T c since 
an excitation with a smaller a is already above the gap 
Ao| sin(2a)|. The relaxation rate 1/r also provides a cut- 
off in the form a m in = 1/tT c . Since the parameter 1/tT c 
is much smaller than all the other cutoff parameters in 
our case, we can write in general 



= max (T/T e ,y/B/Hc 



(14) 



Assuming that U(a < a m i n ) — l/(TaL in ) one can 
insert the truncated potential into Eq. Jll]). Since the 
potential U(a) is concentrated near the nodes one can 
neglect the potential almost in the whole range of a out- 
side a small vicinity of the nodes. Eq. ( O ) has the form 
W = 1+Ae la where the complex constant A can be found 
by matching with the solution in the vicinty of the node, 
a <C 1, i.e., W(a) oc e F ( a \ This provides the boundary 
condition VF(+0) = e 2X W(— 0) across the node at a = 0. 
Here 



2A = F(+0) - F(-0) ~ l/{Ta mm ). 



(15) 



The response function has the same tetragonal sym- 
metry as the underlying system, W(a + n/2) = W{a). 
Together with the above boundary condition, this gives 

1 - e 2A ■ 

W(a) = 1 ^-e la , < a < tt/2, (16) 

1 — ie ZA 

which is to be periodically continued to the rest of angles 
with the period tt/2. We have 



, > 4 1 

(7ff) = l--: 



7T 1 + COtll A 



4 cothA 
IT- (7o) = z ; , -T72T ' ( 17 ) 



7T 1 + coth" A 



Here (• • •} is the average over the azimuthal angle a. 

Magnetic field and temperature dependence. In the low- 
field limit B/H c2 <C (T/T c ) 2 we have A = T c f{TT). At 
not very low temperatures 1/r -C T/T c <C 1, the param- 
eters in Eq. are 



(lo) ~ T C /(TT), 1 - < 7ff ) ~ [r c /(IT)] s 



For T ~ T c , these expressions transform into usual solu- 
tions for a superclean s-wave superconductor 0,0] . For 
very low temperatures T/T c <C 1/r, the parameters be- 
come independent of the scattering time 



( 7o ) = 2/tt, ( 7 h) = 1 - (2/tt). 



(19) 



For fields B/H c2 > (T/T c ) 2 , we have A = 
{l/T)y/H c2 /B. The asymptotics for 1/r 2 < S/i7 c2 are 



{l/T)^H^/B , 1 - ( 1H ) ~ i/ c2 /(r 2 S). (20) 



Eq. (g(j) suggests the magnetic-field dependences of the 
flux-flow conductivity in the form ao oc B~ 3 / 2 . In the 
limit B/H c2 < 1/r 2 we recover Eq. (|l9j). 

Universal superclean limit. For T S> 1 and for low 
temperatures and fields such that 



T/T c and ^ B / H c2 < 1/r, 



(21) 



the longitudinal and transverse conductivities approach 
finite universal values independent of the scattering mean 
free time r: 



(n e + n h )\e\c 
°o = -q (7o) 
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where (70) and (7^} are given by Eq. (|l|). The 
scattering rate affects only the range of temperatures 
and fields, Eq. (21), where this asymptotic behavior 
is established. With lowering the temperature for fixed 
r 3> 1 and B, one expects an increase in ao and, fi- 
nally, a crossover to the regime where both longitudinal 
and Hall conductivities have universal values indepen- 
dent of the mean free path. Together with the con- 
ductivities, the Hall angle reaches its universal value 
determined only by the numbers of electron and holes: 
tan(9# = (tt/2 — l)(n e — rih)/(n e + n^)- The universal 
value of I tan 9h \ is expected to be a local minimum at 
low temperatures due to an increasing ao- 

Experimentally, the sample purity needed for observa- 
tion of the universal behavior of Eq. ( pi"| ) is rather high; it 
requires £/£ » 10 to 100 depending on the ratio E F /T C 
for the particular high-T c superconductor. Nevertheless, 
there are indications that such a regime can be realized 
in practice: the superclean limit was nearly approached 
in Ref. |l3| with T ~ 1 for a presumably d-wave 60 K 
YBaCuO compound |l4| . 

Landau damping. To understand the reason for a fi- 
nite dissipation in the limit of an infinite relaxation time 
t = 00, let us consider the w-dependence of the response 
function W . It is obtained by the substitution ijr — » lu. 
Therefore, A — > —iTroj/4E where 



En 
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2^ 



da 



w (a) 



(23) 



(18) We have 



(7ffM)+j(7oH) = 1- 



2(1 + i) 1 - e-xp(~iTTLu/2E ) 
7r 1 — i exp(— iiruj/2Eo) 



This response function has poles at uik = (4/c + 1)Eq 
where k is an integer. Note that Eq in Eq. (E3h is the 
true interlevel distance in the spectrum of bound states. 
It can be obtained by quantization of the azimuthal mo- 
tion using the Bohr-Zommerfeld rule Q: § Q(a)da — 
2ir{m + 7), where m is an integer and 7 ~ 1. Since 
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Q(a) — —E/uio(a), the quantum-mechanical spectrum is 
E = —(m + j)E Q , assuming that the integral in Eq. J23| ) 
converges due to a finite magnetic field. The poles cor- 
respond to eigenmodes in the quasiparticle distribution 
when the frequency matches the resonance transition be- 
tween the levels with different m. In s-wave superconduc- 
tors, the selection rule for transitions caused by motion of 
an axisymmctric vortex is Am = ±1 with the resonance 
at u> — Eq. In d-wave superconductors other harmonics 
are mixed due to the tetragonal symmetry giving rise to 
resonances at lo = LUk ■ 

When a r „in — * 0, thus Eq — > 0, the poles densly fill the 
u-axis. This leads to a finite density of states and thus 
to a finite Ohmic resistivity at u> = 0: the parameter 



(7o) = -lim Bo ^ Im 

7T 



1 



^ _ g— iiru/2Eo 



1 



ie 



j/2E 



goes to 2 /it as fl — > 0. The attenuation in the zero- 
frequency limit is a special type of the Landau damping 
of vortex motion due to the nodes in the minigap. 

The true interlevel distance Eq introduces a new energy 
scale which does not exist for superconductors without 
gap nodes. For conventional superconductors, there is a 
single characteristic value, the minigap ujq, which marks a 
crossover in the relaxation rate 1/r between two regimes: 
a dissipative viscous flux flow for ujot <C 1, and a dissipa- 
tionless vortex motion for luqt 3> 1. For superconductors 
with gap nodes, there are two critical values: the average 
quasiclassical interlevel distance (luq) and the true inter- 
level distance Eq which can be less than (u>o). Indeed, 
in the superclean limit (uiq)t 3> 1, the true minigap is 
Eq ~ (lu ) a min ; near the nodes, the quasiclassical levels 
approach each other on much shorter distances which can 
only be resolved by quantization of the angular motion. 

The condition of Eq. ((2^) means that E t C 1, A 
finite density of states appears when the true quantum- 
mechanical interlevel distance is smaller than the level 
width. Therefore one has three regimes: (i) Moderately 
clean regime (luq) r« 1 with a high dissipation, (ii) in- 
termediate universal regime Eqt <C 1 with a finite dissi- 
pation at zero u, where the quasiclassical levels are well 
resolved on average but the true minigap is smaller than 
the scattering rate; and (iii) the extreme superclean case 
Eqt 1 when the true levels are well separated and thus 
there is no dissipation. 

The regime with an universal flux-flow conductivity is 
reminiscent of the asymptotic behavior of the current re- 
sponse function in the Meissner state. In the flux-flow 
state, however, the regime of the universal conductiv- 
ity is reached at higher temperatures, T -C T c /T. The 
reason is that the minigap in the spectrum of low-energy 
excitation which are responsible for the vortex dynamics, 
has nodes of a higher order than the gap itself. There- 
fore, the effect of the nodes in the flux-flow regime is 
more pronounced. The universal values of the conduc- 



tivities correspond to an effective scattering rate 1/r of 
the order of the characteristic minigap ujq. This is also 
in parallel with the behavior of the response function in 
the Meissner state, where the effective scattering rate is 
of the order of the characteristic gap Aoo || . 

In conclusion, we have calculated the flux-flow dynam- 
ics in d-wave superconductors at low temperatures. In 
the moderately clean case, A^r/Ep <C 1, the Ohmic and 
Hall components of the flux-flow conductivity are simi- 
lar to those for conventional superconductors. However, 
for superclean systems with A^r/Ep » 1, the magnetic 
field dependence is unusual. Moreover, at low B and 
T, both Ohmic and Hall components approach universal 
values independent of relaxation. A finite dissipation in 
the superclean limit is explained in terms of the Landau 
damping on zero-frequency vortex modes which appear 
due to minigap nodes in the spectrum of bound states 
in the vortex core. Since Aao/Ep is not very small in 
high-T c compounds, the superclean regime is accessible 
in those high purity crystals, where the scaling law for 
the specific heat has been verified fl5| . 
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